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We work out the orthogonality relations for the set of Carniglia-Mandel triple modes which provide 
a set of normal modes for the source- free electromagnetic field in a background consisting of a passive 
dielectric half-space and the vacuum, respectively. Due to the inherent computational complexity of 
the problem, an efficient strategy to accomplish this task is desirable, which is presented in the paper. 
Furthermore, we provide all main steps for the various proofs pertaining to different combinations 
of triple modes in the orthogonality integral. 



I. INTRODUCTION 



The electromagnetic field propagating between boundary surfaces of conducting or dielectric media differs from 
the free electromagnetic field in vacuum in that the boundary conditions imposed by the material media on the field 
alter the structure of the normal modes of the field as well as the state density of electromagnetic field modes. The 
alteration of these quantities becomes manifest when transition rates and energy level shifts of atomic matter exposed 
to the field environment are studied, as these quantities directly enter the expressions of the relevant transition matrix 
elements. 

The complex of problems relating to these ideas is now generally called Cavity Quantum Electrodynamics (CQED). 
One of the principal ideas underlying CQED appeared first in the abstract of paper [Q. CQED as a subject has been 
initiated in with the study of atoms coupled to a single quantized mode of the radiation field in a cavity, followed by 
articles and 0| working on modifications of spontaneous decay rates and radiative level shifts of optical transitions 
in a cavity. Resonators, waveguides, dielectrics, and conducting interfaces modifying atomic radiative properties such 
as decay times, shifts in energy levels, and spectral line shapes have been reviewed in A review of the subject of 
CQED in the presence of metallic boundaries including exhaustive references to original papers was given in Q . 

The mode structure as well as the problem of atomic excitation and photoemission by evanescent waves at a dielectric 
half-space have been derived at first in 0, |^ . There it was found that an appropriate set of normal modes was given 
by so-called triple modes, which are composed of three plane- wave components (incoming, reflected, transmitted). 
The transmitted component of a triple mode can be a travelling wave or an evanescent wave, with exponentially 
decreasing amplitude behind the boundary plane. A variant of these modes has been used in where the radiation 
properties of a quantum harmonic oscillator near a planar dielectric half-space were examined. In , special linear 
combinations of time-reversed and spatially rotated Carniglia-Mandel triple modes have been employed to determine 
spontaneous emission rates and final-state mode densities. A further generalization has been presented in |]ll| , p^ , 
where spontaneous emission in the presence of a dielectric slab bounded by two dielectric half-spaces was studied. 
General quantum optics with the main focus on quantum fluctuations in dielectric media was discussed in ||l3|| . 

One of the central problems in these developments is the determination of an appropriate set of normal modes 
of the free electromagnetic field subject to the given boundary conditions. These normal modes obey the source- 
free Maxwell equations in a dielectric and satisfy the boundary conditions describing either conducting or dielectric 
surfaces. Having solved for a set of modes it is crucial to check their orthogonality, as this property determines 
whether the free field Hamiltonian can be diagonalized in these modes. 

The question of orthogonality of normal modes is simple only in the case of no boundaries, i.e. the field in vacuum. 
As soon as boundaries are imposed, checking for orthogonality can be a tedious procedure. For this reason, actual 
computations proving the orthogonality of modes currently used in the literature are not widely available. For 
example, even in the seemingly simple problem of a dielectric half-space the question of orthogonality of the triple 
modes becomes highly non-trivial. In particular, orthogonality of a pair of triple modes containing a traveling and an 
evanescent component in the same half-space, respectively, is quite counter-intuitive. A brief proof of orthogonality in 
the case of a dielectric half-space has first been given in the appendix of [Q. We feel, however, that there is a demand 
for a more detailed examination of this problem. 
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Our work presented here is intended to fill this gap. We expand on the results presented in the appendix of |j, and 
give a proof of orthogonality for all different cases of triple modes as they appear in the dielectric half-space problem. 
Each proof attempts to highlight the most important intermediate steps as well as the mathematical tools necessary 
to derive them. We do not enter the issue of quantization of these modes, but will end by showing that the free 
field Hamiltonian can be diagonalized in these modes and hence can be written as a sum of independent harmonic 
oscillators pertaining to these triple modes, within a classical context. The completeness of the set of triple modes 
in the space of transverse vector- valued functions has been proven in [ p^ . This, together with our result, shows that 
the triple modes are actually 'good' modes for a canonical quantization of the electromagnetic field in a dielectric 
half-space. 

In this approach, the triple modes appear as traveling normal modes of the free electromagnetic field. The term 'free' 
deserves some explanation here: Strictly speaking, propagation of an electromagnetic plane wave through a material 
medium is the result of the interaction of the fundamental fields E and B with a macroscopic number of microscopic 
sources constituting the bulk matter. In principle, these sources must be incorporated into the dynamics of the total 
system by suitable interaction terms. However, if the frequency band width of the plane waves under consideration 
is far off the resonance frequencies of the bulk matter, any light-matter interaction will be only transient in nature, 
and will certainly not affect the macroscopic state of the bulk matter. In this case, we can treat the bulk matter as a 
passive medium, whose macroscopic properties with respect to the interaction with radiation can be summarized by 
introducing a (space- and time-dependent) refractive index. In the following we will use this approach. 



II. TRIPLE MODES 



The Carniglia-Mandel modes are assumed to propagate in a background which consists of a non-magnetic material 

medium with refractive index n = '^J in the left half-space z < 0, and the vacuum with refractive index n ^ 1 

in the right half-space z > 0, where e, eg are the permittivities of the medium and the vacuum, respectively, and /ip 
is the uniform magnetic permeability which takes the vacuum value everywhere, since the medium is non-magnetic. 
Furthermore, e is assumed to be time- independent. 

Before we state the form of the triple modes, we introduce our notation conventions. Let 63 be the unit normal 
vector onto the boundary plane (i.e. the plane extends in the xy-directions), then any vector V on the boundary 
can be decomposed into components parallel and perpendicular to the boundary plane V — V|| -I- V_l. In particular, 
the wave vector k of a homogeneous plane wave incident on the boundary plane can be decomposed as k =iire3-|-k|| , 
where k • 63 = i^. In the following, this notation will be used extensively. 

The triple modes are obtained by making a separation ansatz for the electric/magnetic fields, factoring out the time 
dependence in the form of e~''^*. A typical electric field will then take the form E (k, x) e^"^*, where the spatial part 
E (k, x) depends on the wave vector k and satisfies an appropriate Helmholtz equation. 

The triple modes can be distinguished by their polarization with respect to the incident plane, which is spanned by 
the incoming wave vector k^ and the normal vector 63. If the incoming electric field E^ is perpendicular to the incident 
plane, the mode is called transverse-electric (TE); this will be indicated by polarization index 1. If the incoming electric 
field lies in the plane of incidence, the mode is called transverse-magnetic (TM), indicated by polarization index 2. 
Each triple mode is composed of an incoming, a reflected and a transmitted plane wave component. We shall label all 
quantities pertaining to these components by subscripts i, r, t, respectively. Furthermore, fields will carry a subscript 
L or R depending on whether the corresponding incoming wave comes from the left, Ki > 0, or from the right, Ki < 0. 

To simplify the subsequent computations we first choose normalization conventions different from those in [Q: we 
assign unit amplitude to incoming electric plane waves in the TE case, and unit amplitude to incoming magnetic 
plane waves in the TM case. Only at the end of our work shall we discuss a different normalization appropriate for 
the diagonalization of the field Hamiltonian. 

The TE modes incoming from the left are labelled by L, the incoming wave vector k^, and the polarization index 
1, and are given by their electric field components 

where the real unit polarization vector e is chosen in the direction of k^ x 83, and the incoming wave is assumed to 
have unit amplitude, k^ is the wave vector of the reflected wave, with Kr — —Ki, and kt is the wave vector of the 
transmitted wave; its component normal to the boundary Kt can be real or imaginary, depending on the angle of 
the incoming wave; in the latter case, the triple mode is called evanescent . The three components of the mode must 
(together with their time dependence e~*"*) satisfy the Maxwell equations on the boundary |l5| as well as a phase 
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matching condition. This determines the reflexion / transmission coefiicients a^, at to be 



For modes incoming from the left, frequency uj and wave vectors k are related by 

kl+Kf =-kl + Kf. = ^loeiu^ , kl + Kf = fioeoLO^ . (3) 

We note that Ki and Kr are always real, but the fact that e > eo allows for Kt to be real or imaginary; in the latter 
case, the transmitted wave is evanescent. The associated TE magnetic field is uniquely determined by Maxwell's 
equations, and is given by 



B^(k.,l,x) = ^ - ^^iJ^^i'x^ ' z>0 ■ 



Similarly, the TE modes normalized upon unit electric field amplitude incoming from the right are labeled by R, 
the incoming wave vector with Ki < 0, and the polarization index 1, and are given by 

V n 1 ^ _ / eate'^'"" , z < , . 

J^^(k„l,xj-|^^^,k..x^„^g.k..x] ^ ^>0 , (5) 

where e has the same direction as k^ x 63. The reflexion / transmission coefficients a^, at have the same form as in 
(^, but the relation between frequency u; and wave vectors is now 

kjj + Kf = k| + K^. = imooJ^ , k| + Kf = ^oeoj"^ ■ (6) 
The associated TE magnetic field is 

. , / at^e"^*"" , z<0 

■t>i?. (Ki, i, X) — < kiXe „ik,.x"i „ k^xe „i;k^.x ' ^ ' 

Now we turn to the TM modes: The magnetic field of the TM modes incoming from the left is 

Bi(k„2,x)- I ^^^^,k..x \ z>0 ' 

where e is a real unit magnetic polarization vector chosen in the direction of k^ x 63. The associated electric field 
follows from Maxwell's equations, and is given by 



kj xe ^ik, «x k^ xe ^ik,.»x z < 

■ ■■■ - i^Moe ' 

Lgikt.x , Z > 



Ei(k„2,x)^ k.xe.g ' ■ (9) 



The relation between frequency and wave vectors is given in eq ) . For TM waves incoming from the right we have 

■D n n \ f ebte'^'"" , z < 

(k^, 2, x) - I ^ ^^,k,.x ^ fo^e.k^.x] ^ ^ > , (10) 



with associated electric field 



{ ktxe ^tkt»x z < 

k,xe^.k"i^r k^g»k..x ' ^ > • (11) 

The relation between wave vectors and frequency is given in eq . 

The reflexion / transmission coefficients br and bt now take a slightly more complicated form in the TM case. 
In order to avoid the need to explicitly distinguish between left- and right-propagating modes, we introduce the 
following notation conventions: Let Hi be the half-space in which the incoming component of the triple mode under 
consideration propagates, and let Si and Ui be the permittivity and refractive index of this half-space; similarly, denote 
Mt for the half-space which contains the transmitted component, with associated permittivity St and refractive index 
rif. 
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Thus, for waves incoming from the left we have Hi = {z < 0} and Hf = {z > 0}, as weh as = e, Ct = eo, 
and rii = Uq, rit = 1; for waves incoming from the right we must interchange i and t in these formulas. With this 
convention, the reflexion / transmission coefficients for the TM waves can be written in a uniform way as 

_ Kje, - Kt/et _^ Kje, 

This result suggests that, in addition to the normal components K — h • n oi the wave vectors, we also introduce 
"reduced normal components" 

Xi = ^ — —K^ , Xt = ^ — —Kt , (13) 
in which case the coefficients (^2|) can be written as 

b — ^ 2.X^i (14) 

Xj^ -\- Xf- Xi -\- x^ 

This notational convention will greatly simplify the computational work in the proofs of orthogonality of the various 
modes. 



III. STRATEGY FOR PROVING ORTHOGONALITY RELATIONS 

We will make use of the following formulas: Let p be a complex number. Then 







dx e^P'' = — , for Imp < ; dx e'P'= = - , for Imp > . (15) 
P J P 

-oo 

As a special case we have 

oo 

lim / dx e^C^-")^ = + t:S (k) ; lim / dx e*(fc+")=^ ^ pi + (k) . (16) 

e^O+O J k e-+0+0 J k 

-oo 

As the computational work required to prove orthogonality of triple modes is considerable, we seek a path of 
minimal effort. It turns out that a good strategy goes as follows: 

We first prove orthogonality of electric fields associated with TE modes, both for co- and counterpropagating 



modes; this is done in section [V. Next we prove orthogonality of magnetic fields associated with TM modes, also for 
CO- and counterpropagating modes; this will be done in section 0. Finally, we show that electric as well as magnetic 
fields associated with one TE and one TM mode always have a vanishing orthogonality integral. This will be done in 
section 0. 



After that, in section VII wc will prove the equations 

J d^x {n^ (x) E* (x) • E' (x) - c^B* (x) . B' (x) } = , (17) 

and 

J d^x {n^ (x) E (x) • E' (x) + c^B {x) • B' {x) } = , (18) 

for all possible TE and TM modes. With the help o f (|l7| ) we can then immediately find the remaining orthogonality 
relations: Applying ( p^ ) to the results from section |IV| proves orthogonality of magnetic fields associated with TE 
modes; while (fjh applied to the results of section ^ yields orthogonality of electric fields associated with TM modes. 
Together with (|1^) we then have completed our task of proving all possible orthogonality relations. 

One might wonder why we choose a strategy as contrived as the one above. The reason is simply that this approach 
minimizes the computational work to be done to prove all the relations discussed above; any other strategy would 
yield the same results, but would require more effort. 
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IV. TE MODES - ORTHOGONALITY OF ELECTRIC FIELDS 
A. Copropagating TE electric fields 

Consider the electric fields of two TE modes 



E(k„l,x) 



eoje 



and 



E' (k^,l,x) 
where both modes are copropagating, i.e. 



sgn {Ki) = sgn {K[) 



X eHt 



X eM, 
X eHt 



Then 



j dPxn^ (x) E(k,;,l,x)* .E' (k;;,l,x) = (27r)^ 5 (k, - k^) 



(19) 



(20) 



(21) 



(22) 



and the result holds for both combinations E^, E'^ and E^, E'^ of subscripts. We see that the square of the refractive 
index serves as a weight function in the integrand. On the right hand side (RHS), rii is the refractive index in the 
common half-space Hi of the two incident components of the triple modes. 
Proof: 

We perform the proof for two modes incident from the left. The case of incident waves coming from the right is 
very similar. If eqs ( p^ , |20| ) are inserted into eq (p2|), the integration over xi and X2 can be performed, which results 
in a delta function 



{2nf 5 {k[ - fci) ,5 (/fc^ - k2) = {2TTf 5(2) (k\^ _ 



(23) 



over the parallel components of the wave vectors. Thus we can replace k|| by k|| in all expressions in the sequel. We 
write this result as 



j d^xn^ (x) E (x)* . E' (x) = (e . e') {2Tif ^f^) (^k|| - k||) x / , 



(24) 



where 



/ dz 



(25) 



and we have omitted reference to the wave vectors k^ and the polarizations 1,2 in E(x). Now we apply theorems 
(|l5Hl6|) to equation (|2^). Since Ki, K[ are always real, the first integral can be treated using theorem (16), resulting 
in two contributions from the principal part and the delta function. The second integral is more involved: If both Kt 
and K[ are real, then theorem ( |l6|) gives 



TT5{K[-Kt) 



(26) 



There are three further cases: K[ real, Kt = iKt; K[ = zKj, Kt real; Kt = int, K[ = in't- Since the incident wave comes 
from the left, the imaginary parts Kj, n't must be positive in each case, and hence the condition that Im {K't — Kt) > 
is always fulfilled, so that we can apply theorem dla): 



dz eHi^'t-Ki) =p 



- K* 



(27) 
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In order to treat all four cases simultaneously, we note that the complex delta function 5 {Kj. — K*) is non-zero only 
in the first case discussed above, i.e. when both and Kf are real. Hence, results ( |2^ ) and ( |2^ ) can be summarized 
into one formula: 

oo 

dze^<^'~^^) =P—±—+nSiK-K:) , (28) 
J^t - 

where it is understood that the (5- function contributes only for real argument, and vanishes otherwise. We can now 
compute expression (eF 



+ n^ta'tatVj^r^ + (39) 
+ nfTT{S{Kl-K,) + a;5{Kl + K,)+a'^d{K^ + K,) + a'XrS{Kl~K,)} + 

where we already have split the principal parts from the delta functions. We now must insert for a^. A long 
computation then shows that all principal parts cancel each other: To see this we use 

{K', - K;) {K[ + A7) = (c.'2 - ..2) ^Q,^ (30) 

to derive 

^2 ^2 

(31) 



{K[ - Ki) {K[ + K,) {K[ - Kt) {K[ + Kt) 

If the last equation is used in ([29|), all contributions from principal parts cancel out. Since both K[ and Ki are positive 
in the present case, the 5 {K[ + ifi)-functions in (|2^) vanish, and wc obtain 

I^nliT{l + a'y,}5{K[-K,)+v?,na'^al5{K[-K;) . (32) 

In this expression, vanishing of the second delta function implies the vanishing of the first. Thus, the only occasion at 
which / is not zero occurs when both Kj. and Kt are real. Then we can express S {K[ — Kt) in terms of 5 {K'- — Ki), 

SiK't-Kt)^4wi^iKl~Ki) , (33) 

where we have used the fact that Kt, Kt as well as Ki, K^ have the same sign. (^2|) now becomes 

/ - 1 1 + a'Xr + a'ta*t ^ | 5 (if^' - X,) . (34) 

Since this expression is non- vanishing only for K'^ = Ki, we can omit the prime's in the curly brackets. Furthermore 
Kt, Ki have the same sign, hence the curly brackets yield 

1 + GrU* + atal — 2 . (35) 
Ki 

Finally, therefore, 

/ = 27^25 (Kl - K,) . (36) 



If ( |36| ) is combined with ( p4| ) we arrive at the result (|22^ 

The same result holds for two copropagating modes incident from the right: Again, the square of the index of 
refraction (x) emerges as weight in the integral, while now on the RHS nf is equal to 1. 



B. Counterpropagating TE electric fields 



In a similar way we prove the orthogonality of two counter-propagating TE modes. It is sufficient to assume that 
the field E is incoming from the left, while E' is incident from the right: 
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Consider the electric fields of two counter-propagating TE modes 



and 



Ek(k'„l,x) 



z < 
z>0 



where K, > and K' < 0. Then 



Proof: 



J d^x (x) Ei (k„l,x)*.E'^(k'„l,x) = 



Following the same course as in the previous section we find that 

j d^xn^x) E(x)* .E'(x) = (27r)^ (e . e') .5(2) (k,, -k||) x / 
where the quantity / now takes the form 



(37) 



(38) 



(39) 



(40) 



I ^ Hq J dz 



j dz 



(41) 



Since the result is non- vanishing only for k'| = k||, we have e • e' = 1. Integration over z and splitting into principal 
part and delta functions using cqs (^-^6|) gives for the principal-part-contribution to ( ^l|) 

T _ AiKi 1 



^princ 



{K[+K[){K,+K;)li? 



-Kf 



^1 K'J^-K-^ ^ "-2 K'-^-Kf^ 



(42) 



while the delta contributions to ( pl| ) all vanish, since Ki and K[ have now opposing sign. Thus, eq (^) can be 
rewritten as 



AiK,Kl{K', + Kt) 



Using the fact that 



and 



*2 _ ^2 (-^12 7^2 



- K. 



£l 



2 _ ^2^2 



nn = — n 



(43) 

(44) 
(45) 



the curly brackets are seen to vanish, which proves the result. 



V. TM MODES - ORTHOGONALITY OF MAGNETIC FIELDS 

Here we prove analogous relations for TM modes. Since we start with magnetic rather than electric fields, we 
normalize the incoming magnetic field in such a way that the magnitude of the associated electric field E^ is 1, 



|B.| = M|E,| = i 

OJ c 



Hence, the incoming magnetic plane wave takes the form 

1 



c 



, |e| = l 



(46) 



(47) 
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A. Copropagating TM magnetic fields 



Consider the magnetic fields of two TM modes 

B(k„2,x) = 

and 

B' (k^,2,x) = 
where both modes are copropagating, i.e. 



r. ^ 



•X 



X (Mt 



X ee, 

X (Mt 



sgn (K,) = sgn {K[) 



Then 



/ d3xB(k„2,x)*.B'(k'„2,x) = — (2^)' (5 (k, - kO , 



(48) 



(49) 



and the result holds for both combinations B^, B^ and B/j, B^j of subscripts. 

Remark : It is important to note that, in contrast to the electric case (|2^), in the magnetic case no refractive 
index appears in the integral (and neither on the RHS of the above equation) . The same is true for counterpropagating 
modes. 

Proof: 

We proceed similar to the proofs given above. In a first step, we perform the integration over the coordinates 
Xj| = {xi,X2) in (^), giving the intermediate result 



j d^x B (x)* . B' (x) = (e . e') S'^^^ (k,, - kj,) x / 



(50) 



Due to the delta factor we have e • e' = 1. As before, the expression / can be written as the sum of two parts 

^delta + I] 

vanishes as in the previous computations, /princ = 0. The delta contribution is 



delta -i- ^princj involving the delta functions and the principal values, respectively. The principal part contribution 



/delta = t: {b;b',6 {K[ ^Kt) + S {K[ - K,) + b'X^ {K[ - Ki)} 



(51) 



where the first term involving 5 {K[ — Kt) is present only if both Kt and K[ are real. If cither of them is imaginary 
then according to the previous discussion the first term vanishes. 
Next we use the fact that the modes are copropagating, 



to obtain 



sgn {Ki) = sgn {Kl) = sgn (Kt) = sgn (Kf) 



(52) 



(53) 



provided that both Kt and K[ are real. The quantities X are defined in eq (p^). 

Now, three cases must be distinguished: (1), both Kt and K[ are real; (2), both Kt and K[ are imaginary; (3), one 
of them is real and the other one is imaginary. We first deal with case (3) : We have the equation 



Kl^'-^\K?~KJ) . 



In case (3), the LHS is always nonvanishing, and so is the RHS. But this means that certainly 



(54) 



(55) 



hence all delta functions involving Ki and K[ vanish in (|5^). Furthermore, the first delta function involving Kt and 
K[ vanishes naturally in this case, as was explained in section [V. It follows that the orthogonality integral must 
vanish in this case, and so eq (Efl) is confirmed. 
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Now we turn to case (1), where Kt and are real. Here we find on using (|53|) that 



/delta = 7r6 {K[ - i^,) J + 1 + 



(56) 



Direct computation now shows that the curly brackets yield 2, and therefore 

/delta = 2^<5 {K'^ - K,) . (57) 

(|^) together with (|^) now confirms (^9|). 

Finally, we consider the case where both Kt and K[ are imaginary, and hence are associated with evanescent modes. 



Now we have = —Kt, and (51) gives 

/delta = 7r(S(/f,'-/^0{l + |6r|'} • (58) 

But 

|b.r= '^''^*'' =l , (59) 

due to the fact that Xi is real, and Xt is purely imaginary, X^ = —Xt- Thus, again we find (|5^) for /delta, which 
confirms (^9|). ■ 

B. Counterpropagating TM magnetic fields 

This case proceeds exactly analogous to the previous results, so we quote the result without explicit proof: The 
magnetic fields of two counterpropagating TM modes, 

sgn{Ki)^~.sgn{K';) , (60) 

are orthogonal, 

y"d3a;BL(k„2,x)*.B^(k:,2,x) = . (61) 

VI. ONE TE- , ONE TM MODE 

Finally, we show that a pair of modes such that one member is TE and the other one is TM, is always orthogonal. 
This holds for both electric and magnetic fields: 

J d^x n^{x)'E{k,,l,x)* •E' {k'„2,x) ^ (62) 

y"d3a;B(k„l,x)* .B'(k^,2,x) = . (63) 

These results hold for all combinations of subscripts LL, LR, RL and RR. 
Proof: 



From the results above, we know that each of the integrals J in eqs (|62|, |63|) has the form 

J = (2^)' (5(2) (kii -kJi) X / , (64) 



which follows from performing the integration in (|6^) over the two coordinates X|| . Because of the delta factor, both 
wave vectors lie in the same plane of incidence. First, consider the electric case in (^), and assume that E is TE, 
and E' is TM. Then E is perpendicular to the plane of incidence, while E' lies in the plane, hence the LHS of ( |6^ ) 
vanishes. ~ The same argument immediately applies to the magnetic case. I 
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VII. A THEOREM LINKING ELECTRIC AND MAGNETIC ORTHOGONALITY 

Now we turn to prove the following equations, 

J d^x{n^ (x) E* (x) • E' (x) - c^B* {x) •B'{x)}=0 , (65) 

J d^x{n'^{x)'E{x)»'E'{x)+c^B{x)»B'{x)} =0 , (66) 

valid for all combinations LL, LR, RL and RR of subscripts as well as for both TE and TM polarizations. These 
equations will complete our discussion of orthogonality of triple modes. 

Remark: Note that the first of these equations involves one complex conjugate of an electric and a magnetic 
field, respectively, while the second equation contains no complex conjugate fields. 

Proof: 



We first prove (|65 
We start by introducing the explicit time dependence of the modes. 



E = E (a;) e"^'^* , B = B (a;) e"*"* 
E' = E' (x) e-"^'* , B' = B' (x) e 



(67) 



These time-dependent fields must obey Maxwell's equations, from which we derive 



V • [E* X B'] = i l^'^E* • E' - ujB* • B'l 



(68) 



This equation contains a time dependence expi {to — uj')t on both sides. Integrating over t and uj' yields 

-V.[E*(a;)xB'(x)]= J^E*(x).E'(x)-B*(a;).B'(a;)l , (69) 



where E (a;) etc are again triple mode functions. In order to prove ( |65[ ) we therefore must show that 

j d^x\/» [E* (x) X B' (a;)] = , (70) 

for any two triple modes E (x) and B' (x). To see this we first note that the spatial dependence of the integrand can 
be written as 

E*(a;) xB'(x) =C(z)e'^"*(''"-''ii) , (71) 

where C (z) is given by the expression 

C (z) = e-^^ii'K"*'") [E* (x) X B' (a;)] . (72) 

A brief investigation shows that C(z) indeed depends only on z, as suggested by the notation. Insertion into ( [70| ) 
yields 

J d^xV* [E* (x) X B' (x)] = (27r)^ (J^^) _ j^^^^ j ^q^q (^^j (73) 
It follows that we must show the vanishing of the surface integral 

j dz[d,C{z)]=Q , for kj|=k|| . (74) 



We see immediately from eq (|7^) that, if E is TE and B' is TM, then C and hence the integral (|73| ) vanish as required. 
To show that the same holds if E is TM and B' is TE, wc take the complex conjugate of eq (|69|), 



-V . [E* (x) X B' (.x)] = |-i-V . [E'* (x) X B (x)]l 
UJ I ^ J 



(75) 



Since by assumption, E' and B are parallel, the RHS of ([75|) vanishes, which proves the above statement. 
Thus, we are left to prove the following four cases: 
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1. E and B' TE, copropagating; 

2. E and B' TE, counter-propagating; 

3. E and B' TM, copropagating; 

4. E and B' TM, counter-propagating. 

Due to the S^^^^ (k'^^ — k||^ factor in (|73|), we can assume that both wave vectors lie in the same plane of incidence. 

This imphes, in particular, that in cases (1) and (2), the magnetic field B' lies in the common plane of incidence, 
while E is perpendicular to it; and in cases (3) and (4), the electric field E lies in the plane of incidence, while B' is 
perpendicular. 

Let us now treat case (1): For the sake of simplicity, we assume that both modes propagate from the left to the 
right, i.e. Ki, K[ > 0. We decompose the divergence in eq (173) into parallel and normal derivatives. 



V.(E* X B') = V|| • (E* X B')|| +a3(E* X B')3 . (76) 

A computation shows that integration of the first term on the RHS of (|7^) over the spatial coordinates vanishes 
because of property ([75|). In the second term, all terms involving delta functions vanish. Evaluating the principal 
value contributions gives 



/ d^x u'ds (E* X B')3 = (27r)' |^ _ k,, j lu' ■ {K[ {1 - < -f- < - «} - ^^^^ 



On account of 

K 



^r^r = -^a^a^ , (78) 



we find that the RHS of (77) must vanish. The statement that integral ( [7 3) ) vanishes is therefore confirmed, for case 
(!)■ 

Now we briefly discuss the remaining three cases: Case (2) proceeds along the same lines as case (1). To prove 
cases (3) and (4) we start with 

V« (E* X B') = -V« (B' X E*) . (79) 
On the RHS, we express E by the associated magnetic field B, which is TM, 

„ k; X Bi „ kr X Br „ kt X Bf , , 

E, = , E, = , Et = . 80 

Lo^iQCi uj^ioei Lo^jiQtt 

Subsequent evaluation of the integral then yields a factor {2tt)^ 5^"^^ ^ky — k||^ times an expression 

■^f {Xi,Xt,X'^,X^,br,bt,b'j.,bf), where / as a function of the arguments X and b has the same functional form 
as the analogous expression in the TE case as a function of respective TE quantities K and a, 

/ {X,,Xt, Xi Xi br, bt, 6;, b[) ^ f {K,, Kt.Kl K'„ a„ at, a[) . (81) 

Now, the functional form of the quantities b in terms of X is the same as the functional form of a in terms of K, as 
can be seen by comparing eqs @i and (p^. It follows that we may transfer all conclusions derived previously for the 
TE case to the TM case by making the replacements indicated in the last formula (^l|). In particular, vanishing / on 
the RHS of ( pTl ) implies vanishing / on the LHS. In this way, we have traced back cases (3,4) to cases (1,2), which 
finishes the proof of (|6^). 

— Now we can turn to the proof of eq (^6|). This proceeds analogous to the previous equation (|65|), hence we give 
only an outline: By using the explicit time dependence ( |67[ ) of the modes we can derive the relation 

- V • [E (x) X B' (x)] ^ i ( ^E (x) • E' (x) -I- B (a;) • B' (x) \ (82) 



UJ - ■ - - - , ^2 



between triple mode functions, which is analogous to eq (|69|). To prove ( |66[ ) we therefore must show that 

j SxV • [E {x) X B' {x)] = , (83) 
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which is analogous to ([70|). The left hand side of the last equation can again be cast into a form similar to (|73|), 

Jd^xV*[E*ix)xB'{x)] = {2TTf6'-^^(k\+k\\^Jdz[d,C{z)] , (84) 

except that the sign of the second wave vector in the argument of the delta function now differs. It follows that we 
must show the vanishing of the surface integral 

J dz [d,C (z)] = , for kj| = -k|| , (85) 



which is analogous to (74). The remainder of the proof proceeds along the same lines as the proof for (65). This 
finishes our proof of eqs (|65|, |66| ) . 



The remaining orthogonality relations now follow from (|65|) : Suppose that E and E' are electric fields associated 
with TM modes; then orthogonality follows from orthogonality of the magnetic fields as proven in section [v| and (|6^). 
On the other hand, suppose that B and B' are magnetic fields associated with TE modes; then their orthogonality 
follows from orthogonality of the associated electric fields as proven in section |^and again eq (^5|). This finishes our 
proof of orthogonality relations. 

VIII. NORMALIZED MODES 

In the developments above, we have normalized the electric field of TE modes upon unit amplitude of the incoming 
plane wave, and the magnetic field of TM modes upon a value numerically equal to 1/c for the incoming plane 
wave. This choice was motivated by simplification of the subsequent computational work. Now we finally introduce 
a normalization which is best suited for the purpose of diagonalization of the free field Hamiltonian, as performed 
below. The new normalization will affect only TE modes; in this case, eq ( p2[ ) motivates that we divide these modes 
by the refractive index Ui of the medium in the incoming half-space H^. If the TE mode under consideration is 
incoming from the right, the associated refractive index is actually 1, so that these modes are unaffected. Thus, the 
new normalization effectively changes only the TE modes incoming from the left. 

The new TE modes are thus defined as 

Ei(k„l,x)^^ Jk..x ^ \>Q ' (86) 



-e ■ at ■ e 

no ^ 



and 



with 



K, - Kt 2Ki 
at 



K, + Kt ' K, + Kt 

Furthermore, B (k^, l,x) is defined to be the magnetic field associated with (|8^, ^), where each of the plane wave 
components of B is given in terms of the plane wave components of E according to 

„ ki X Ei „ kr X Er „ kt X Ef , , 

B^=- , B,-^ , Bt^- . 89 

LO Ll) LO 

Formulas ( ^ , ^ ) are valid for both left- and right-incoming modes. This defines the normalized TE modes . 
The normalized TM modes are given by 



B, fk- 2x1-/ ^^[^^'"'" + ^^'--^'''^"] ' 
±iL(k„AxJ- j ie . . e'kt.x ^ ^>o 



(90) 



where the reflection and transmission coefficients follow from eqs (|12|-|1 

K,, ~ nlKt , 2K,, 



Ki + n„Kt Ki + ngKt 
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and for modes incoming from the right, 



B« (k 2 X) - ^ K'"" + ^''^*^] ^ '-^ (92) 



with 



- if* , 2K, 



The plane wave components of the associated TM electric fields are given in terms of plane wave components of the 
TM magnetic fields: For waves incoming from the left we have 

„ ki X Bi „ kr X Br „ okf X Bt , , 

El« = 2^ ' Ei, = ^ , Eii = -c2^ i , (94) 

while for waves incoming from the right, 

_ X B, „2 kr xB, kt x Bt 

-tiflj — ~C , i^Hr = -C , Ji^m = 2 ■ \^^) 

All orthogonality relations discussed before can now summarized: 

J d^x (x) E* (k, s, x) • E' (k', s' , x) = {2nf Sss'S{k' - k) , (96) 

J d'x W (k, s, x) . B' (k', s', x) = {2TTf ^Sss'S (k, - k^) , (97) 
for both combinations LL and RR, and 

j (fx {x) (k, s, x) • E'^ (k', s', x) = , (98) 

y"d3a:B2(k,s,x).B'fl(k',s',x) = . (99) 



Furthermore, eqs (tea, 66) continue to hold. 



IX. EXPANSION OF ARBITRARY SOURCE-FREE FIELDS IN TERMS OF TRIPLE MODES 

The set of triple modes is complete with respect to the source-free fields, as has been shown in This means 
that every source- free electric / magnetic field can be uniquely expressed as a linear combination of triple modes. For 
an arbitrary free real electric field E (x, t) we can make the ansatz 

E (x,t) = (100) 
= I J d3fc^£(fc) [^^^(k)Ei(k, s,x) (k)E2(k,s,x) e^' 



I J d^kJ2 [uRs (k) Efl (k, s, x) - (k) E|j (k, s, x) e'"*] 



where 



£:(fc)^J^iM_ . (101) 

V 2eo {2irf 
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The associated magnetic field is determined by Maxwell's equation V x E = — B, and has the expansion 

B (x,t) = (102) 

= i J d^kY, ^ (fc) [uLs (k) Bl (k, s, x) e-^"* - ul, (k) B^ (k, s, x) e*"*] 

K<0 

2 

+ I / d3fc^f(fc) [7.j^,(k)Bfl(k,s,x)e-'"*-uJj,(k)BJj (k,s,x)e^"*] . 
1 



We note that the amplitude factor £{k) is the same as in formula (100) 



X. EXPANSION OF THE HAMILTONIAN IN TERMS OF TRIPLE MODES 

The Hamiltonian of the source-free radiation field is based on the classical electromagnetic field energy 



H=^ Jd'^x [n^ (x) (x,i) + c^B^ (x, t)] , (103) 



with real fields E and B. If we insert expansions (IOC, 102) into (103), we obtain the 
mode expansion of the Hamiltonian 

2 

H= [ d'kJ2^ l^Ls (k) ul, (k) + ul^ (k) (k)] 

s— 1 
2 

+ / d^fcE^Ks(k)«L(k) + «L(k)"i^.(k)] . (104) 



K>0 



To arrive at this expression, we had to employ all orthogonality relations (96-391) and (pq, |66|) 



In our current discussion, the coefficients ulr,s are c-numbers, as we are still in the classical domain. In deriving 



the mode expansion (104), we have been careful, however, to retain the order in which terms appear, so that (104) 
can be immediately transferred into the quantum regime, by interpreting uls (k), (k) etc. as annihilation and 
creation operators of triple modes. 

XI. SUMMARY 

We have worked out the orthogonality relations for the set of Carniglia-Mandel triple modes, which provide a set 
of normal modes for the source-free electromagnetic field in a background consisting of a dielectric half-space and the 
vacuum, respectively. The inherent computational complexity in this kind of problem creates a demand for efficient 
strategies to accomplish this task. In this paper we have provided this strategy, and furthermore, we have presented 
a comprehensive discussion of the proofs for the various distinct cases in the orthogonality integral in order to serve 
as a compilation of technical details for related problems in the field of optics with possibly even greater technical 
complexity. 
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